arXiv: 1507.06488v2 [hep-th] 30Jul2015 


arXiv:1507.06488 


Entanglement Temperature With Gauss-Bonnet Term 


Shesansu Sekhar Pal“ Sudhakar Panda^ iU 

“Department of Physics, Utkal University, Bhubaneswar 751004, India 
^Institute of Physics, Sachivalaya Marg, Bhubaneswar 751005, India 

shesansu lATIg maifcom pandaHliopb.res.in 
Abstract 

We compute the entanglement temperature using the first law-like of thermodynamics, 
AE = TentASEE, up to Gauss-Bonnet term in the Jacobson-Myers entropy functional in any 
arbitrary spacetime dimension. The computation is done when the entangling region is the 
geometry of a slab. We also show that such a Gauss-Bonnet term, which becomes a total 
derivative, when the co-dimension two hypersurface is four dimensional, does not contribute 
to the hnite term in the entanglement entropy. We observe that the Weyl-squared term does 
not contribute to the entanglement entropy. It is important to note that the calculations are 
performed when the entangling region is very small and the energy is calculated using the 
normal Hamiltonian. 
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1 Introduction and summary 


The recent study of entanglement entropy has drawn a lot of attention because of its remark¬ 
able similarity with the black hole entropy [1]. The entanglement entropy is dehned as von 
Neumann entropy: S = —Tracey {paLoq pa), where pA is the the reduced density matrix. 
For a spherical surface of radius, TZ, the reduced desnity matrix is dehned by tracing out 
the degrees of freedom that sits inside this radius, which is the complement of A. In which 
case, the von Neumann entropy depends on the area of the spherical region, S ~ TZ^. This 
von Neumann entropy is interpreted as the entropy seen by an observer sitting outside the 
radius 7Z. 

In the context of gauge-gravity duality [2], a prescription, strictly speaking a conjecture, 
is suggested by Ryu and Takayanagi (RT) [3] to calculate the von Neumann entropy in 
the gravitational side. This von Neumann entropy is called as entanglement entropy. The 
RT prescription suggests to consider a co-dimension two spatial hypersurface in a way such 
that its boundary coincides with the boundary of the region that we are interested in and 
then hnd out the area of the minimal surface. Finally, the entanglement entropy is the 
ratio between the area and where G^i is the Newton’s constant, (for recent reviews, 

see HI El E])- A claim of the proof of the RT conjecture is provided in different spacetime 
dimensions and with different entangling region in [H E] . The subadditivity nature of the 
entanglement entropy is shown in [9]. 

In this paper, we shall report the result of a holographic calculation of the entanglement 
entropy of a given region upon inclusion of the terms up to Gauss-Bonnet term in the 
Jacobson-Myers (JM) entropy functional [TU] in any arbitrary speactime dimension. The 
shape of the entangling region that we are interested in is that of the slab type. To recall, 
the prescription that we shall follow to carry out such a calculation is that given by Ryu and 
Takayanagi (RT) [3], of course without the higher derivative term. 

The JM entanglement entropy functional up to four derivative term is 


AGmSee — d' 
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det{gab) 1 + ^iR{g) + A [R {g) — 4:Rab{g)R°' {g) + Rabcd{g)R 


yabcd 


( 1 ) 

where Ai and A are unknown real coupling constants and are dimension full. A derivation of 
the JM functional starting from the Einstein-Hilbert action with the higher derivative terms 
are given in [II11I21E3]. For our purpose, we do not require the full derivation of it. 

The induced metric on the co-dimension two hypersurface is denoted as gab = daX^dbX^C 
where Gmn is the d + 1 dimensional bulk spacetime geometry and the hypersurface is de¬ 
scribed by X^. The precise form of the hypersurface that follows (HUS] 
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which is essentially the equation of motion associated to the held and 
whose precise form /Cf^ = dadbX^ — 'y^f^dcX^ + daX^dbX^Tf^j^, where 7 ^ 5 , and hf^jy 
the affine connections dehned using the induced metric gab and the bulk geometry, Gmn, 
respectively. This particular form of the hypersurface holds good irrespective of the shape 
and size of the entangling region. 

Let us note that without the higher derivative terms the equation of the hypersurface 
is derived earlier in [16] and is called as the extremal surface. In what follows, we shall be 
interested in the strip type entangling region only. 

In this paper, we shall compute the correction to the expression of the RT entanglement 
entropy by considering such higher derivative terms in the JM functional. This essentially 
means we are computing the entanglement entropy in the hnite ‘t Hooft limiij^ In fact, we 
shall be doing the calculation to the linear order in the couplings Ai and A, for simplicity. 
The result reads for the bulk geometry as AdS spacetime with radiu^ Rq as 
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where the entangling region is taken as slab type. The slab that we are interested in is (d — 1) 
dimensional. Along one direction, it is extended from 0 to size i and along the rest of the 
directions from —L/2 to L/2. It is more properly defined in eq([^. The UV regulator is 
defined as e. 

As found by RT and [18], the entanglement entropy has a divergent piece, which comes 
from UV and a finite piece from IR. However, there are few salient features that are worth 
mentioning: 


^It simply follows from the AdS/CFT correspondence [2]. The quantities Ai and A are dimension full 
and the only allowed quantity that come is the size of the AdS spacetime, R. By the above mentioned 
duality R ^ where is and A are the string length and the ’tHooft coupling, respectively. It means 

Ai ~ and A ~ 


^The size Rq is the solution to the bulk equation of motion with higher derivative term whereas R is 
without m- The relationship is as follows: Rq = where foo is the positive real root of the following 

v/~ 

cubic equation 1 — /oo + _ {d-2){d-3Xd-4){d-5)A what follows, we shall express 

everything in terms of Rq only, for simplicity. 
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a) This result of the entanglement entropy makes sense only for d > 3. For d = 3, it is 
easy to notice that the hnite piece does not receive any correction at hnite ‘t Hooft coupling. 
However, the divergent piece can receive corrections. 

b) The power of the UV regulator, e, is independent of the value of the ‘t Hooft coupling, 
i.e., it remains the same for inhnite as well as hnite value of the ‘t Hooft coupling. However, 
the coefficient of it depends very well on the value of the ‘t Hooft coupling. 

c) The entanglement entropy depends on the quantity i via power law type and this 
behavior is not changed even in the hnite value of the ’t Hooft coupling. However, i does 
depend of the value of the ‘t Hooft coupling through the turning point r* to be followed 
latter. 


d) The power of L, the size of the slab along the other spatial directions, remains same 
both in the divergent as well as in the hnite piece. This is because, it does not play any role 
in the determination of the extremal surface and comes as an over all factor in the integration 
of the entanglement entropy. This is simply because of the translational invariance along 
these directions. 

e) Note, for d = 5, the Gauss-Bonnet term becomes topological, which means it won’t 

contribute in the determination of the hypersurface, This fact is rehected in eq(12), 

to be followed. Moreover, such a term becomes a total derivative and does not contribute to 
the hnite term in the entanglement entropy. 

Doing a small perturbation around the AdS spacetime geometry described by a parame¬ 
ter, m, is shown to respect the hrst law-like of thermodynamics AE = T^ntASEE [IH], where 
AE corresponds to the energy of the excited state and ASee is the entanglement entropy. 
The entangling temperature depends on the size of the entangling region £. More impor¬ 
tantly, this is shown in the limit of small entangling region, m£^ <C 1 as in [12]. Some other 
interesting studies on the hrst law-like thermodynamics is done in [20] EU |22l |23l EU ES], EE] . 

In this paper, we shall, test the hrst law-like of thermodynamics, which reads as AE = 
TentASEE, for the low-lying excited states with higher derivative corrections in the entangle¬ 
ment entropy functional for any arbitrary spacetime dimensions. For the slab type entangling 
region, the entangling temperature takes the similar form as before: T^nt = cji, where c is 
a constant and depends on the couplings. With the higher derivative term as written in the 
JM functional 


c = 


2(^2-!)/ 2{d-2){d-3) 2(d-2)(d-3)(d-4)(d-5)(3d2-6d + 2) \ 
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For d = 3, the quantity c is universal in the sense that it is a pure number and does not 
receive any corrections in the hnite ’tHooft couplings. 
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The paper is organized as follows: In section 2, we give the compntational details of the 
entanglement entropy with higher derivative terms, np to Gauss-Bonnet term. In section 3, 
we show the hrst law-like of thermodynamics with higher derivative terms and compute the 
entanglement temperature. In section 4, we re-visited the concavity and the specific heat 
and then conclude in section 5. 


2 Entanglement Entropy 

Let us consider the form of the entanglement entropy functional as written in eq([^ for the 
strip type entangling region. In this case the extremal hypersurface that follows is given in 
eq@ with 

/Cf, = dad,X^ - + daX^dtX^Tl,^, (5) 

where 7 ^^ and Tf^jy the affine connections defined using the induced metric gab and the 
bulk geometry, Gmn, respectively. The indices a b, c etc run over the codimension two 
hypersurface, whereas M, N, S etc run over the entire spacetime. 

2.1 Slab/strip type entangling region 

Let us consider the following background geometry possessing the translational symmetry 
along the temporal and spatial directions along with the rotational symmetry, with diagonal 
form as 


= GMNdx^dx^ = -gtt{r)dt^ gxx{r){dxl H-h grr{r)dr‘^. ( 6 ) 

The induced geometry that follows on a codimension two hypersurface is 

UjT 

X* = 0 , = x“, X" = r(xi), r' = — 

dx\ 


ds'^_i = gabda'^da^ = 


gxxij'') T Qrrij'^ 


' dr 
dxi 


dx\ + gxx{r){dxl ^ - Vdxl_^). (7) 


In order to carry out the explicit computations, we shall use the following form of the 
slab type entangling region 


0 < xi < £, -L/2 < (x 2 , • • •, Xd-i) < L/2. 


( 8 ) 


Various coordinate invariant quantities involving the Riemann tensor on the codimension 
two hypersurface take the following from 


R = 


(d-2) 


^9lx{9xx + 


2r'^gxx9'xx9'rr -id- ly^gxxg'lx - ^r"glJxx - 


5 






rCab-^ — 


RabrdR 


abed 


4'r QxxQrrQ XX 5)r QrrQ 


J2 

XX 


id-2) 


id - 2) {grrg'iy^ + gxxi2gly‘^ + g'rrg'xxr"^ - 2grrr'^g"xx) - 


^QgLiSxx + grrr'^Y 
2‘g‘lxi'^'‘^ 9 " XX + g'xxr")^ + (^d — 4:)grrg'lx'i"'^ + ‘ 2 gl^ir'‘^g"xx + g'xxr") + 


gxx^ 


n 


id 5)g XX 9 rr9 xxf + 2grrf' g 
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Using these informations, we can compute G-B = RabcdR^’’^'^'^ — 4i?“^i?ab + and comes 


?ab i 


as 
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{d - 2){d - 3){d - 4) 
^QgLidxx + grrR^y 


/2 /2 


(d - + gxxr'^(id - 13 )c/'L - 4 g'rrg' xxt'^ + 


r'V'.. +8^/i(r'V'.. + d'.y') 
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In passing we must mention that upon doing the calculation of the Weyl-squared term 
dehned as: WeyR = RabcdR^’’^^ — ^j^R'^'^Rab + (^d- 2 )(d- 3 ) vanishes identically for the strip 
type of entangling region. Hence, to conclude the Weyl-squared term does not contribute 
anything to the entanglement entropy functional. 

Substituting all these into the equations of motion eq([^ and denoting x[ = ^ gives 

^ ( 9 %^x'i _ {d-2){d-3)gj~x'lg'l^ ^^ ^ jd - 2){d - 3){d - 4)(d - h)gRf g'tx^'iA 

dr 1 g^r gxxX'l 4:igrr gxxRl)^R 16((7rr + 

( 11 ) 

Solving the equation to leading order in the couplings 


x\{r) = 
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where c is the constant of integration. It is chosen to take c = gxx (r*), because we have 
opted the following boundary condition at r = r*: the quantity Limitr^r*a:'i = 00 , diverges. 
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Substituting it into the entanglement entropy functional gives 


2GnSee = L 
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Let us evaluate the above integral by consider a specific example. For our purpose, we 
shall take the example of AdS spacetime. For this purpose, we put the boundary to be at 
r = 0 ^ 

dsl_^_^ = ^ + dxl H-h (14) 


The result reads as 
2GnSee = 




^ (d-l)(d-2)^ , (d-l)id-2)id-3){d-A) ^ 

(d - 2V-^ d Rl + Ri - 


^L-‘-^Rt'^(wR)) 6 , (d-2)(d-3) , (d-l)(d-2)(d-3)(d-4)(d-5) ,\ 

(d - 2 )ri-^ r I fig ‘ ( 2 d - l)fi3 )’ 

(15) 


where the integral over r is performed from the the UV cutoff, r = e, to IR r = r*. Now we 
can use the following relationship between £ and r* to leading order in the couplings 

(d-1)(d-S), , (d-i)(d-1)(d-S)(d-i)(d-b) 

' r(^) b'- fig - (2^-i)fig 

(16) 


This follows from eq(12), upon doing the integration. On substituting this into the entan¬ 
glement entropy as written in eq(15) gives us eqg. 


2.2 For (i = 5 

It is interesting to note that for a specihc dimension, d, the Gauss-Bonnet term does not 
contribute anything to the equation of motion. This happens when d = 5, in fact, in this case, 
the Gauss-Bonnet term in the entanglement entropy functional becomes a pure topological 
term. In fact, it can be expressed as a total derivative. The purpose of this subsection is 
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to show that this total derivative term does not contribute anything to the hnite piece but 
does to the divergent piece of the entanglement entropy. 

For d = 5, the Gauss-Bonnet term can be expressed as 
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f^dxr^det{gab) {R\g) - + Rabcd{9)R"'’"\9)) = -\r x 

^^ 9xx - ‘^9xx9Di9rr + 9xxX'^) + 9xx9'xxi9'rr + 9'xxXf + 2g^^x'^x'[)) 
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Using the bulk as AdS spacetime structure and evaluating it gives the desired result 
as written above. So to conclude, we hnd that inclusion of a total derivative term in the 
entanglement entropy functional does not contribute anything to the hnite term of the en¬ 
tanglement entropy. 


2.3 Scaling symmetry 

Let us demand that under the following scale transformation of the coordinates, the metric 
components transformations as 

x^ —)■ Xx^, Gmn A ‘^Gmn ( 18 ) 

for which the length between two points in the bulk geometry remains invariant. This 
property holds for the AdS spacetime apart from other symmetries. On the codimension 
two hypersurface, we also want 

(t“ -)■ A(t“, gab X~‘^gah, gabdcr'^da’' -)■ X°gabda°-da''. (19) 

Then it simply follows that the Lorentz scalar quantities made out of the Ricci curvatures 
behaves as 

R{g) ^ ^ R^Ua)R°‘"‘‘(a) ^ >~'‘R<^(g)R‘‘^‘'‘(.9)- 

( 20 ) 

It is very easy to check that the JM entanglement entropy functional eq(|3 does not 
scale under the above transformation. It is because of this scaling symmetry the size of the 
entangling region comes as an over all factor in the computation of the entanglement entropy 
functional. One can notice that each term in the following expression 

1 + X^Rig) + A [R\g) - ARab{g)R^\g) + RabcMR^^’^^'^id)) (21) 
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does not scale as said above. In fact, we expect it not to scale for a scale invariant bulk 
geometry. It is also expected that the size L should not enter in the computation of this 
expression. Hence, it is quite natural to expect that the size of the entangling region should 
come as an over all factor. Note that both L, ^ scale as (L, £) —)■ A(L, t). This means the 


turning point r* also scale linearly in A by virtue of eq(16). It also means that the quantity 
c scales as c —)■ 

From our studies of thermodynamics, it is well known that if we scale the volume V —>■ AV 
and energy E —)■ AE, then the entropy scales as S{V, E) —)■ S{AV,AE) = AS{V, E). This 
is just the homogeneity of order one property obeyed by the entropy. 

Note for the AdS geometry the energy vanishes which means the only relevant function 
is S{V). Recall, for the slab type entangling region, the volume V = In order to 

have the desired scaling for the volume, we must scale {£, L) —)■ A^(f', L). This means 
A = A^. In which case, the hnite part of the entanglement entropy does not scale but the 
same cannot be said about the singular part. As it is a bit ambiguous. 


3 Small Fluctuations 


In this section, we shall check the hrst law-like of thermodynamics by considering small 
fluctuations of the bulk geometry along the lines of [in] but now with the higher derivative 
terms. In order to check such a law, the entanglement entropy will play the role of the 
entropy and the change in energy will be the energy of the low lying excited states and then 
their ratio will give us the analogue of the temperature which will be called as entanglement 
temperature. 

For the following type of fluctuation of the AdS geometry 
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where m is the fluctuating parameters. We can do the necessary calculations very easily 
using the expressions for the entanglement entropy as given in eq([l^. In what follows, we 
shall assume that <C 1 [I9]. This condition can as well be re-written as 1, so 

as to keep terms to linear order in the parameter, m. Now, to satisfy such a requirement 
means we must take r* very close to UV because the boundary is at r = 0. 


Substituting these metric fluctuations into eq(13) and carrying out the necessary r inte¬ 
grals from UV to IR results in 
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Now, we want to express r* in terms of i to leading order in m and it comes as 
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Now, we can use such a relation to re-express the entanglement entropy in terms of i as 
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There follows the expression of the change in entanglement entropy: ASee = See{'^, £) — 
SEE{'m = 0, £), which is hnite and does not depend on the UV cutoff, e. 

Energy: The energy of such a perturbed geometry can be calculated using the expression 
of the energy momentum tensor as given in 1271 with the prescription given in |28j . 

AM = J d‘^~^x\Jdet{aij) N u^u^Tmn, (26) 
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Tmn = 
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where K is the trace of the extrinsic curvature. The other quantities are defined as 


Jmn = \{2KKmlK^ + KlsK^^Kmn-2KmlK^^Ksn-K^Kmn), 
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The quantity, i?, acts as a regulator and we can expand R = Rq + XRi + KR 2 to linear order 
in the couplings. The sizes Ri and R 2 will be determined by demanding that Tu becomes 
finite as we approach the boundary. Or in the limit of m —?• 0, the Tu component should 
vanish as well [2H]- The quantities Ri ^2 are 




(d-2)(d-3)(d-4)(d-5), 


(29) 


where Rq is related to R as Rq = R/ vToo and foo obeys the following equation: 1 — /oo + 
(d-2Kd-3) _ (d-2)(d-3Kd-4)(d-5) ^y^ ^ g shall uot express any of our 

results in terms of R, for simplicity. 


which in our case using N = u* = 1/y/^ and Jdet{a)ij = gxi , gives 


AM = J d'^-^x g. 


^ Tt 
2 ^ 
XX 
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(30) 


The expression of the various components of the extrinsic curvatures are Ku = ~ 2 ^^) Kxx 

/K- and K = „ ^'*1_ + Substituting all these along with the AdS spacetime into 

the above integral results in 
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(31) 

The bulk couplings A, A are related to the couplings appear in the entanglement entropy 
functional are A = ^ and A = A im. It is easy to see that for vanishing m, the energy 
vanishes, which is true for the AdS spacetime [28]. Hence AE = E{m). 
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Having got all the desired expression for the change in the entanglement entropy and the 
change in energy means, we can evalnate the entanglement temperatnre, T^nt = ‘^E/ASee = 
j. Upon doing the calculations, we obtained eqQ, which is true for any d > 3 when the 
entangling region is of the slab type. One of the interesting point is that the quantity c is 
independent of the parameter m, even though the energy and the entanglement entropy do 
depend on it. 

There exists a difference in the way the entanglement temperature is calculated in com¬ 
parison to that in the black hole physics, even though there exists a first law-like relation of 
thermodynamics in both cases. Note that the entanglement temperature is not the inverse 
periodicity associated to any shrinking one-cycle as in black hole physics, because there does 
not exist any on the co-dimension two hypersurface. In fact, it takes the following form 




= dxl 


tfxx 


+ 


{d-2){d-3)g'lAc^-gt-^) 


\ 2c grrdxx 

{d -2){d- 3){d - 4)(d - 5)g%{gi-^ - c^) 


Ai + 


2\2 




A I -l- gxx{dx\ dx‘^_i),{32) 


where Note that for the AdS type background geometry the co-dimension 

two geometry is not conformal to flat space for any dimension d > 3. Even though, it is 
not conformal to flat space but it respects the scaling symmetry as discussed previously for 
m = 0. 


4 Concavity & Specific heat 

Concavity: It is suggested that the entanglement entropy obeys the strong sub-additivity 
property and this followed by looking at the concavity of entanglement entropy (TH]. It 
means the concavity automatically implies the sub-additivity of entanglement entropy. The 
property concavity suggests the second derivative of the entanglement entropy with respect 
to the size of the system should be negative: ^ < 0. 

In our case, the entanglement entropy at hnite ’tHooft coupling takes the following form 


as in eq(25). 


See = Sf 


a 


T d—2 j->d—l I T d—2 jyd—l f}2 

L/ rtn ~\ -— u JXn 


"'0 


(33) 


where is the UV regulator dependent part and the value of a, [3 can be obtained by 


comparing with eq(25), which are essentially positive. Demanding the concavity condition 
gives 

EE ^ n , n ^ ~ ~ 2) 


d'^S, 


dP 


< 0 




2mP 


(34) 
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Remember, we work in a limit for which 1, and this means for hnite a this condition 

is fulhlled automatically. Hence, according to [19] the entanglement entropy at hnite ’tHooft 
coupling obeys the strong sub-additivity property. 


Specific heat: It is suggested in [12] that the specihc heat that follows from the calcu¬ 
lation of the entanglement entropy becomes positiv^ For unit dynamical exponent, the 
temperature goes as T^nt = c/£. In which case, the specihc heat dehned as 


d{/\SEE) 

^ ent 


dT 


= -2 


ent 


m(3 cVRl 


d-l 


(35) 


ent 


where V = iL'^ ^ is the volume of entangling region under study and ASee is the last term 


in eq(33). 


There are a couple of comments in order: 

(a) The volume is a function of temperature, V = ^ and it goes as inverse of the 

entanglement temperature. Form which it follows trivially that = — 1. 


V dV 


(b)The negative form of the heat capacity (because /3 > 0) implies the system under 
study is unstable, which is in agreement with the fact that we are trying to understand 
the properties of an excited state. Just to compare, the temperature dependence of the heat 


capacity for a black hole with spatial horizon in dj-1 dimensional AdS spacetime, C ~ 
where Th is the Hawking temperature. 

(c) In thermodynamics, upon scaling the volume V —>■ AV and energy E —)■ AE^ makes 
the entropy to scale as S' —)■ AS. This is due to the homogeneity condition of the entropy. 
However, such a simple relation does not hold in the case of entanglement entropy. The 

d 

change in entanglement entropy, ASee goes as ASee —t A'^ ASee- 

Generically, the holding of the second law of thermodynamics, ASee A 0 , suggests the 
positivity of the specihc heat. However, in the present case, it looks like that the excited 
states does not obey the second law of thermodynamics. 


5 Conclusion 

To conclude, we have calculated the contribution of the higher derivative terms up to Gauss- 
Bonnet term in the JM entanglement entropy functional by considering the background 
spacetime as AdS. To carry out the calculations we have considered the entangling region is 
of the slab type. At the end, the form of the hnite piece of the entanglement entropy up to a 

L^d—2^d—1^2 — d 

sign goes as SEEifinite) n-u -^- —/(Ai, A), where the functional form of the couplings 

/(Ai, A) is very difficult to predict and follows by doing the calculations. The dependence 
on L, Ro and i follows from translational invariance and dimensional analysis. However, the 

feel that such a calculation should be revisited. 
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absence of the information related to the entangling region like L and £ in /(Ai, A) follows 
simply from the scale invariance. Moreover, we observed that the Weyl squared term does 
not contribute anything to the entanglement entropy functional. 

For d = 5, i.e., when the bulk spacetime is 5 + 1 dimensional, it is easy to see that the 
Gauss-Bonnet term in the JM entropy functional becomes a pure topological term. In fact, 
this piece becomes a total derivative term. Moreover, upon calculating we hnd that it does 
not contribute to the hnite term but does contribute to the UV regulator dependent term 
in the entanglement entropy. 

Recently, it was suggested in [2H| by adding a different kind of exotic term to the en¬ 
tanglement entropy functional, which is a total derivative term in the bulk, and was shown 
that such a term does contribute to the hnite part of the entanglement entropy. So, the 
future question of interest would be to hnd under what condition does a total derivative 
term contributes to the entanglement entropy, generically? 

After doing a small perturbation described by a parameter, m, around the AdS spae- 
time, we hnd that the change in entanglement entropy up to an overall couplings goes as 

The energy associated to such a 


d—2 pci—1 /p2 




Cjv 
rriL' 


d—2 ryd—l 


^— Now, if we we demand that 


ASEE{'m,i) = SEE{'m,i) - SEEim = 0,£) 

geometry up to an overall couplings goes as AE - ^ 

there exists a hrst law-like of thermodynamics as in m then the entanglement tempera¬ 
ture dehned as Tent = = f) where the quantity c depends on the couplings and the 

dimension of the spacetime. 

In this paper, we have worked in the limit when the entangling region is very small, 
mi'^ -C 1, and checked the hrst law-like of thermodynamics using the normal Hamiltonian 
for the slab type geometry. However, in [32] and [33] , the hrst law-like of thermodynamics is 
checked using the modular Hamiltonian for the ball-shaped entangling region. It is certainly 
very interesting to check this for other type of entangling geometries, which we leave for 
future studies. 
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